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Abstract 

The thermodynamic properties of S ^ 1 ferromagnetic chains with an easy-axis single-ion 
anisotropy are investigated at arbitrary temperatures by both a Green-function approach, based 
on a decoupling of three-spin operator products, and by exact diagonalizations of chains with up 
to ./V = 12 sites using periodic boundary conditions. A good agreement between the results of both 
approaches is found. For the 5=1 chain, the temperature dependence of the specific heat reveals 
two maxima, if the ratio of the anisotropy energy D and the exchange energy J exceeds a charac- 
teristic value, D/J > 7.4, and only one maximum for D/J < 7.4. This is in contrast to previous 
exact diagonalization data for comparably small chains (N ^ 7) using open boundary conditions. 
Comparing the theory with experiments on di-bromo Ni complexes the fit to the specific heat yields 
concrete values for D and J which are used to make predictions for the temperature dependences 
of the spin-wave spectrum, the correlation length, and the transverse magnetic susceptibility. 

PACS numbers: 75.10.Jm; 75.40. Cx; 75.40.Gb 
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I. INTRODUCTION 



In low-dimensional spin systems^ the interplay of quantum and thermal fluctuations and 
the effects of spin anisotropics on the thermodynamics are of basic interest. Whereas for 
Heisenberg antiferromagnets quantum fluctuations occur already at T = 0, in ferromagnets, 
possibly with an easy-axis anisotropy, quantum fluctuations exist at nonzero temperatures 
only. The study of systems with ferromagnetic exchange couplings, e.g., of the quasi-one- 
dimensional (ID) S = 1 ferromagnet CSN1F3 with an easy-plane single-ion anisotropy— is 
also motivated by the progress in the synthesis of new low-dimensional materials, such as the 
frustrated S = ~ quasi-lD cuprates.^ Likewise, the magnetic behavior of LaMn03, a parent 
compound of the colossal magnetoresistance manganites,— may be described by an effec- 
tive spin S — 2 model with a ferromagnetic intraplane and an antiferromagnetic interplane 
coupling, where neutron-scattering experiments 5 yield evidence for a pronounced ferromag- 
netic short-range order (SRO) in the paramagnetic phase and for an easy-axis single-ion 
anisotropy. To provide a good analytical description of SRO and of the thermodynamics at 
arbitrary temperatures, the standard spin-wave approaches cannot be adopted. Recently, 
the Green-function equation of motion decoupling of second order and the Green-function 
projection method with a two-operator basis, respectively— i 7 ? 8 ! 9 ! 10 ! 11 ! 12 ! 13 have been success- 
fully applied to quantum spin systems with anisotropies of different kind, where most of 
the previous work was devoted to S = ~ systems and only in Refs. [D and |9| low-dimensional 
isotropic 5^1 ferromagnets have been considered. However, the study of anisotropic S ^ 1 

nrin] 

magnets is of current interest . 14 i 15 i 16 i 17 In Refs. Il5lll6lll7l the S = 1 ferromagnet with an easy- 
axis single-ion anisotropy was investigated in the random phase approximation (RPA) for 
the exchange term. By this approach the paramagnetic phase and its SRO properties cannot 
be described. Therefore, a second-order Green-function theory of SRO for anisotropic S ^ 1 
models, going one step beyond the RPA, should be developed. 

As a first step in this direction, we present a theory for the 5^1 ferromagnetic chain 
with an easy-axis single-ion anisotropy described by the model 

H = - J -Y J S % S ] -D"£{Stf (1) 

denote nearest-neighbor (NN) sites] with J > 0, D > 0, and Sf = S(S + 1). The 
choice D > is motivated by our emphasis on the temperature dependence of the specific 
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heat, which reveals a double maximum being more pronounced for D > than for D < 0,— 
and by the comparison with the experiments on di-bromo Ni complexes^ which may be 
described by the model (JTJ) with D > 0. 

Furthermore, we perform exact finite-lattice diagonalizations (ED) of S — 1 chains with 
up to iV = 12 sites using periodic boundary conditions which are critically analyzed in 
relation to the ED results by Blote^ for the specific heat using open boundary conditions. 

The paper is organized as follows. In Sec. II the second-order Green-function theor 
for th e model (JTJ) is developed extending previous approaches for 5 = ~ 
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12 13) to S ^ 1 and rotation-invariant methods for S ^ 1 (Refs 
an easy-axis on-site anisotropy. To test the theory, in Sec. Ill the limiting cases J = 
and D = are considered in comparison with exact and ED results. The effects of spin 
anisotropy are explored in Sec. IV. The spin-wave spectrum, the spin susceptibility, and the 
specific heat are investigated with the focus on the condition for the existence of two maxima 
in the temperature dependence of the specific heat. Moreover, the theory is compared 
with available experimental data, and predictions for some relevant quantities are made. A 
summary of our work can be found in Sec. V. 

II. GREEN-FUNCTION THEORY 

The dynamic spin susceptibilities Xq^i^) = ~((Sqi S- q ))u (yn = H — , zz; — 7r ^ q ^ ir), 
defined in terms of two-time retarded commutator Green functions,— are determined by 
the projection method . 10 ' 11 ' 12 Taking into account the breaking of rotational symmetry by 
the single-ion spin anisotropy we choose, as in Refs. and Q, the two-operator basis 



A u = (Sq,iSq) T [y = +,z). Because the model considered has up-down symmetry with 
respect to S? — > — Sf, we have (Sf) = 0. Neglecting the self-energy, the matrix Green 
function ((A; A + )) w — [u — M / M _1 ] _1 M with the moment matrices M = ([A, A + ]) and 
M' = ([iA,A+]) yields 

= - 2 lw (^) 2 = M^/M^, (2) 

where = ([iS^ , and Mq 3 ^^ = ([— , — iS1 q ]). The first spectral moments are 

given by the exact expressions 

M+~ = 4Jd(l - cosg) + 2D[3C* Z - S{S + 1)], (3) 
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M zz = 2JC+-(l-cosq). (4) 

The correlation functions C n = \C+~ + C z n z and C v / = (S%S%) = j^T C v ^e iqn with 
C£ M = (SgStg) are calculated by the spectral theorem,— analogous to Ref. |l2|, as 

M UfM 

c 7 = ^ l + 2n H^ + D 7i ( 5 ) 

D7 = ^((S»;S» q ))P, (6) 

where n(c<j) = (e w / T — and ((•••;*' denotes the anticommutator Green function. 
The on-site correlators Cq^ are related by the sum rule 

Ct + C zz = S(S + 1) (7) 

which follows from the operator identity Sf = S^S^ — S z + (S z ) 2 and (S z ) = 0. 

To obtain the spectra uf^ in Eq. (j2J) in terms of two-spin correlation functions we approx- 
imate the time evolution of the spin operators ~S q in the spirit of the schemes proposed 
in Refs. O^lsOni 11 12 . That is, taking the site representation the products of three spin 
operators in —S\ are expressed in terms of one spin operator. Then, the projection method 
neglecting the self-energy becomes equivalent to the equation of motion decoupling in second 
order. 

In — we decouple the operators along NN sequences as^ 

S+S}Sr = a+-(S+S r )S+ + a+-(S+SnSf. (8) 

Here, following the investigation of the isotropic ferromagnet,^ the dependence on the 
relative site positions of the vertex parameters (cf. Ref. Q) is neglected. 

For S ^ 1, in —Sf there appear products of three spin operators with two coinciding 
sites which we decouple as proposed in Refs. HandEl 

SfSjS} = (SrS+)S+ + \ + -(S+Sr)S+. (9) 

Furthermore, for contains the term D 2 Ai with 

Ai = S+(S Z ) 2 + 2S z SfS z + {SlfSf. (10) 

For S = ~ we have Ai = 0, and for S — 1 we get Ai = (Ref. 13) using the relation 
= S'St (Ref. Q. To obtain a r ea S0 nab,e approbation of A, for S > 1, we 



calculate exactly the average {A i S i )(T) at T = and T — > oo. We obtain (A i S i )(0) = 
(2S - l) 2 (5+5r)(0) with (S+Sr)(0) = C+~(0) = S and lim T _ >00 (A5'-) = ±[4S(S + 1) - 
3] liniT- >0 o(5'j f ST") wr th h m T^oo(S'j + S'r) = ^(S 1 + 1). Due to those results, for S > 1 we 
approximately replace Ai by 

A, = V (2S - 1) 2 S+ (11) 

with rj(T = 0) = 1 and r/(T -> oo) = ^j0jw - Note that Eq. (JUJ) holds exactly for 
S = 1 with r/(T) = 1. Considering the ratio R = t)/Cq z , forS 1 = 2 (3) we have lim^oo -R = 
0.23 (0.09) as compared with R(Q) = S~ 2 = 0.25 (0.11). Accordingly, for 1 < S ^ 3, i?(T) 
depends only weakly on temperature. Neglecting this dependence, i.e., taking R(T) = 
R(0), T](T) in Eq. (jllj) may be calculated in a reasonable approximation as 

V(T) = ^QTCO- (12) 



In — S zz we adopt the decouplings (cf. Refs. 
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and [J 
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SfSjSr = a zz (S+S-)S z , (13) 
SrS*S+ = \ zz (S+Sr)S z . (14) 



Finally, we obtain the spectra 



[1 - cos q) { A + - + 4J 2 a + - d{l - cosg)} + (cu+~) 2 , (15) 



A+- = J 2 {S(S +1) + C ZZ + 2A + -C 1 + 2a + -(C 2 -3C 1 )} 

+ 2DJ{2X + -C ZZ + 3C ZZ - S{S + 1)}, (16) 

(u+-) 2 = 2DJ{S{S + 1) - 3C ZZ + X + -(2C ZZ - Of")} + rj{2S - 1) 2 D 2 , (17) 

{u zz ) 2 = (1 - cosg){A 2Z + AJ 2 a zz C+-{\ - cosg)}, (18) 

A zz = 2J 2 {S(S + 1) - C zz + a 2Z (C+- - 3C+-)} + 2J(J - 2D)A Z2 C+-. (19) 

To calculate the correlation functions from Eq. (jSJ), in particular the term given 
in Eq. (jBJ), we follow the reasonings of our previous paper.— We obtain D+~ = 0, because 
WgS Q 7^ for D 7^ 0, and D zz = ^ n Cfi 9i o. Then, the longitudinal correlation functions are 
calculated as 



with C zz = jj J2n^n an d C zz given by the first term in Eq. (jSJ). Note that the term C zz 
describes long-range order in the infinite system. In previous wor k 8 ' 10 ' 11 ' 13 such terms are 
introduced by hand and interpreted as condensation parts. 

By Eqs. . . (fTHJ) . and (fTTJJ) we get the longitudinal static susceptibility 

X\ z = Xo Z {l + ±a zz C+-(/\ zz y\l - cosg)}- 1 , (21) 

where Xq z = j&zz ■ Expanding the denominator for small q up to 0(q) we obtain the 
correlation length £ zz = \/2a zz C^ / A zz . 

The uniform static susceptibilities xX ma y be also expressed in terms of C^. From 
Eqs. (J2J and (jSJ) with lim^ C^ = Tlim q ^ x z q z , we get 

Xo=^J2 C n> ( 22 ) 

n 

which agrees with the general formula of thermodynamics in the case (Sf) = 0. For Xo~ we 
obtain 

xt-=gJ2 C n~ (23) 

with g = 2{u£~[l + 2n(u£~)]}- 1 , following from Eqs. ©, ©, ©, and (fTTHl. 

For large temperatures, T 3> the static susceptibilities and the structure factors 

C q ^ are related by 



At very high temperatures, in Eq. (f2"2|) only the n = term may be taken into account, and, 
with C +- = 2C^ = lS(S + 1), we get the Curie law = 2 X ZZ = 2S(S + 1)/3T. 

To provide a better comparison of the Green-function theory with ED data, it is useful to 
consider the theory also for finite systems with periodic boundary conditions. For a ring with 
an even number N of spins we have the discrete q values = with — y rij ^ y — ^ 
and y + 1 correlators with ^ n ^ y. In the calculation of C zz according to Eq. (j2TiJ) 
we must take care of the term C zz which, for finite N, is finite at arbitrary temperatures 
and is given by 

1 / N/2 ~ l \ 

C zz = \c zz + C z N z /2 + 2 J2 C n\- (25) 



n=l 
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FIG. 1: J 



limit: Specific heat for S > 1. The Green-function theory for S = | (dashed) 



and 5 = 2 (solid) is compared with the exact results for S = | (dotted) and S = 2 (dot-dashed), 
respectively. 

Then, it turns out that two equations for C^ z are linearly dependent. As additional equation, 
we use the expression of Xq z m Eq. ()21j) in terms of according to Eq. ()22|). i.e. 



2C+ 
JA 1 



(26) 



III. LIMITING CASES 



To test the approximations made for S ^ 1 in addition to those for S — ^, in particular 
the decouplings © and (JHJ), where X U ^(S = \) = 0, and the replacement (fTTj) with Eq. (fT^j). 
we first consider the limiting cases J = and D = 0. 

In the J = limit, by Eqs. ©,©,©, and (02D we obtain 

2S{S + 1) 



3 + v ^(25-l)[l + 2n( 



(27) 



where Wg" = ^frj{2S — 1)D and is calculated by Eq. (fT2*|) . The longitudinal on-site 



correlator Cg z is obtained from Eq. (J7J). At T = and T — > oo we get Cq z = S 2 and 
^(^ + 1), respectively, agreeing with the exact results. Figure [T] shows the specific 
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heat for S 



| and S 



' C (S - 
and T£ 



1.4D for S 



2 derived from Cq 2 , where the temperatures of the maximum 
ic 

2, respectively. This yields a justification for the 



T%(S = |) = 0.97D and T^(S = 2) = 1.6D nearly agree with the exact values Tg = 0.83D 



| and S 



approximations (fTTj) and (JT2J). The result for 5 = 1, not depicted in Fig. [T] for clarity, 
shows qualitatively the same temperature dependence as that for S > 1; it is exact, because 
Eq. (|11|) becomes the exact relation = S^. 

In the D = limit, we have C+~ = 2C* Z , a Ufl = a, X v " = A, A^ = A, and = u q . 
The vertex parameter a(T) is determined by Eq. (J7J), Cq z = 5(5 + l)/3. To derive an 
equation for A(T), we first consider the long-range ordered ground state with £ _1 (0) = 
corresponding, by Eq. (J2U, to A(0) = 0. Then, by Eq. (JT5J) we have u q = 2J ^/2aC{ z {\ — 



_____ 

cosg) and, by Eq. (|20|) . C* z = \ -£-8 n ,a + C zz . Taking into account the exact result 



CS (0) = |5 2 we get a(0) = § and A(0) = 2 — | (cf. Ref. tJ). At non-zero temperatures 
there is no long-range order, i.e. C zz = and A > 0. To improve the approximation of 
Ref . 0, A(T) = A(0), we first derive the exact high-temperature series expansion of C\ z up 
to 0(T- 2 ), 

CZi 
1, 



1ZZ 

,ex 



s(s + di V J u 2 



a j yf-lh) + °^- < 28) 

Expanding Eq. (p0)l for n = 1 and n = up to (^(T" 1 ) and using, for n = 0, Eq. (|25j) we 
obtain 

"5(5 + 1) 



Cf 



«o^, (29) 



<T = {l " T^ 3 + AS ( S + !)( A o " »o)]^} , (30) 

where «o and Ao are the lowest orders in the expansions of ot(T) and A(T), respectively. The 
comparison with the exact results Eq. (J2%|) and Cq z = 5(5 + l)/3 yields 

ao^l, A = 1-3(45(5 + 1)]" 1 . (31) 

The result c_o — 1 confirms the general suggestion (cf. Refs. that the vertex 

parameters a approach unity at high temperatures. Considering the ratio Q = A/a, for 
5=1,2 and 3 we have lim^^oo Q = 0.63, 0.88, and 0.94, respectively, as compared with 
Q(0) = 0.67,1, and 1.1. Accordingly, for 1 ^ 5 ^ 3, Q(T) is only weakly temperature 
dependent. Setting Q(T) = Q(0), A(T) may be calculated in a rather good approximation 
as 

A(T) = ^2-4W). (32) 



3 V 

In Fig. 121 our results for 5 = 1 are plotted, where a remarkably good agreement of the 
Green-function theory for iV = 12 with the ED data is found. This justifies the decouplings 
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FIG. 2: D = limit: Nearest-neighbor correlation function Cj 1 , longitudinal spin susceptibility 
Xq z (upper inset), and specific heat C (lower inset) for S = 1, as obtained by the Green- function 
theory in the thermodynamic limit (solid) and for a finite system with N = 12 (dashed) in compar- 
ison with the ED data for N = 12 (o). The Green-function results for obtained for N = 12 
and N — > oo agree within the accuracy of drawing. 

(JUJ) and (fT3j) with A calculated by Eq. (p?2*|) . Considering the specific heat, the temperature of 
the maximum in the theory, TIr = 0.9 J for both N = 12 and iV — > oo, only slightly deviates 



221 agrees 



from the ED result = 0.8 J. Note that in the semiclassical approach of Ref. 
with the ED value. Concerning the uniform static susceptibility in the thermodynamic limit, 
we have (see also Ref.|9j) lim T ^ Xo z T 2 / J = ^S 4 . This low-temperature behavior, \l z °^ T~ 2 , 
qualitatively agrees with the result of the renormalization-group approach of Ref. Q, but 
quantitatively deviates from the finding liniT^o XqT 2 / J = 1.58S 12 . 23 



24 



for the 



Finally, let us compare our theory with the Green-function approach of Ref. 
5 = 1 antiferromagnetic Heisenberg chain. There, instead of the decoupling (|14j). the 



left-hand side is rewritten as 5,- S z Sf 



Q— C+ Qz 
°i °j °j 



S^S^. The first term is decoupled 
analogous to Eq. (fblj) as a(S^ 5^ + )5| , whereas the second term yields a contribution to i Sf. 
This results in a gap A in u* z at q = which, for D = 0, is given by A = J. In Ref. |2J, 
A is interpreted as a Haldane gap. As we have verified, A is independent of 5. However, 



for 5 



2' 



for example, there is no Haldane gap. That means, the gap A is an artefact 



of the approach of Ref. |2J employing commutation before decoupling. According to our 
experience (see, e.g., the Green-function theory for the t — J model 2 ^) such a procedure 



should be avoided. Furthermore, we argue that the approach of Ref. |2J yields A ^ for 
the S — 1 antiferromagnet also in higher dimensions and for the S ^ 1 ferromagnetic chain. 
Concluding, contrary to the reasonings of Ref. 13, the Haldane physics cannot be captured 
by the second-order Green-function theory. 



IV. EFFECTS OF SPIN ANISOTROPY 

To complete our Green-function scheme for the model with D > and J > (here- 
after, we set J = 1), the four parameters a Ufl and X utl have to be determined. In the ground 
state, for D > we have the exact results 

C^(0) = S5 nfi ; C^ z (0) = S 2 (33) 

so that CWO) = S 2 . By Eq. © we get C+~(0) = ^ V ^e^ n and, comparing with 

E q . dnnj), 

M ? +- = 25w+". (34) 

Inserting M+~ and cu^ given by Eqs. (JHJ) and (fTo")) to (|T7j) with 77(0) = 1 [see Eq. (JTTJ)] and 

comparing the coefficients in Eq. (j34j) in front of (1 — cosg) n (n = 2 and or 1), we obtain 

a+-(0) = l; A+-(0) = 1-^. (35) 

Considering finite temperatures and suggesting limr-^oo (T) = 1 (see Sec. Ill), we put 
a + ~(T) = 1 because of a + ~(0) = 1. Following the reasonings in the D = limit, for the 
ratio Q + ~ = A + ~/a + ~ we assume Q + ~{T) = Q + ~(0), i.e. A + ~(T) = A + ~(0). The parameter 
a zz (T) is calculated from the sum rule (J7j). Concerning the remaining parameter X zz and 
the ratio Q zz = X zz /a zz , it turns out that Q zz has very different values in the T — > 
and T 00 limits. Therefore, we adjust A ZZ (T) to the ED data for C ZZ (T) which are 
depicted, for S = 1, in the inset of Fig. El Thus, we have a closed system of equations for 
seven quantities (Cq^, C^, a zz , \ zz ) to be determined self-consistent ly as functions of 
temperature. 

As a first test of our approach, in Fig. |3] the NN correlation function for S — 1 is 
plotted, where a very good agreement with the ED results is found. The correlator Cq~ (not 
shown) also agrees very well with the ED data. As can be seen, we have < 2C\ Z ] that 
is, due to the easy-axis anisotropy the transverse correlations are suppressed as compared 
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FIG. 3: Nearest-neighbor transverse spin correlation functions C{ for S = 1 at D = 0.1, 1,5, 10, 
and 25, from top to bottom, showing the Green-function (solid) and ED results (o, N = 12). The 
inset shows the ED data for the nearest-neighbor longitudinal correlation functions Cf 2 at the same 
values of D with D increasing from bottom to top, which are used as input to the Green-function 
theory. 

with the longitudinal correlations. The maximum in the temperature dependence of 
indicates the crossover from Ising-like to Heisenberg-like behavior, where the maximum 
position increases with increasing D. 



A. Spin waves 



At T = 0, by Eq. (|34j) with Eqs. (j2J) and (|33j) we obtain the spin- wave spectrum 

uj+~(0) = 2S(l-cosq) + (2S -1)D (36) 

with the spin-wave gap Uq~(Q) = (2S — \)D. Let us point out that the dispersion (j3T?j) 
agrees with the result obtained by the RPA and the Anderson-Callen decoupling (see, e.g., 



Ref. ll6J) given by B { = S^S? + SfS? = 2(S Z ) {l - ^ ■ [S(S + 1) - Cf]}5+; putting, at 
T = 0, (S z ) = S and CgJ = S 2 so that B t = 2S - 1, the spectrum ® results. In the 
RPA approach of Ref. [Hi where the D term for S = 1 is treated exactly, we calculate 
Xq~(^) = -2(u - ujq)' 1 (correcting a misprint in Eq. (51) of Ref. 17) with uj q given by 
Eq. flUl) with 5=1. 

Let us compare Eq. ()36|) with previous spin-wave theories. The generalized spin-wave 

11 
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FIG. 4: Spin- wave spectrum for S = 1 at T = 0, 1, and 5, at q = it from top to bottom. The 
inset shows the spin-wave gap lOq at D = 0.1, 1, and 5, at T = 1 from bottom to top. 

theory by Becked for 5=1, which extends the Holstein-Primakoff transformation to two 
sets of Bose operators treating the single-ion anisotropy exactly, yields cJq~~(0) = D, in 
agreement with Eq. (|36|). Contrary, in the ordinary spin- wave theory (with only one Bose 
operator a*), (S z ) 2 with S z = S — rii and rij = a+cij is approximated as (S*) 2 = S 2 — 2Srii 
neglecting the nf term. This yields the wrong result lo^~ = 2SD violating the condition 
Uq~{S — |) = 0. Note that such an approach was used to fit the inelastic neutron-scattering 
data on LaMnOa on the basis of an effective spin model with easy-axis single-ion anisotropy.^ 
From our results we conclude that this fit should be reconsidered by means of an improved 
theory. 

In Fig.|3]the temperature dependence of the spin- wave spectrum for S = 1 is shown, where 
a spin correlation-induced flattening of the shape with increasing temperature is observed. 
The spin-wave gap as function of temperature exhibits a minimum and approaches the high- 
temperature limit lirnT->oo^o"Cr) = ^o"~(0). In the paraphase (T > 0) with SRO, well- 
defined spin waves exist, if their wavelength is much smaller than the correlation length, 
i.e., if q 3> (^ zz ) _1 - To estimate the validity region of the spin-wave picture, in Fig. El the 
inverse correlation length is plotted. For D = we get lim^o £ ZZ T = S 2 (cf. Ref. |9|) which 
nearly agrees with the result of the renormalization-group approach,— lim-r^o C, ZZ T = 1.145* 2 . 
For D > the low-temperature behavior of £ zz is quite different. By Eq. (121)) we have 
^ zz = y / 2a zz I 'A zz and {i zz )~ 2 Xo — /a zz with a zz = a zz C^, where the numerical 
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FIG. 5: Inverse correlation length for S = 1 at D = 0.1, 1, 5, 10, and 25, from top to bottom (solid), 
compared with the D = limit (dashed). 

evaluation yields a finite value of d zz as T — > 0. Because (xo 2 ) _1 (0) = 0, (£ z ^)~ 2 approaches 
zero as T — > much stronger than Cf~ and {Xo)~ l (compare Fig. |5] with Figs. H3 and 
EJ). Correspondingly, the easy-axis anisotropy drives the paraphase at low temperatures 
close to long-range order. As can be seen from Fig. the validity region of the spin-wave 
picture, q ^> (£ 2Z ) _1 , shrinks with increasing temperature, where predominantly high-energy 
magnons may be observed. 

B. Spin susceptibility 

The spin anisotropy results in a qualitatively different temperature dependence of the 
uniform static susceptibilities Xo~ anc ^ Xq z i as can be seen from Fig. |U1 Note that the ED 
calculation of Xo~ requires a small magnetic field in the x direction so that the S — 1 data 
can be obtained only for N = 8. 

The transverse susceptibility xt ~ (Fig-EK) reveals a maximum at T*, where increases 
with D (right inset), in very good agreement with the ED results. For small anisotropies 
(left inset) a pronounced finite-size effect is observed, where the theory for N = 8 agrees 
well with the ED data. The temperature dependence of xt~ ma J be explained as follows. 
The anisotropy-induced longitudinal SRO (cf. Fig. results in a spin stiffness against the 
orientation of the transverse spin components along an external field perpendicular to the 
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FIG. 6: Transverse (a) and longitudinal (b) uniform static spin susceptibility for S = 1 at D = 
1, 5, 10, and 25, from left to right, obtained by the ED (o) for N = 8 (a) and N = 12 (b) and by the 
Green-function theory (solid). In the insets on the left-hand side, for D = 0.1 the Green- function 
results for N — > oo (solid) and N = 8 (dashed) are compared with the ED data for N = 8. In the 
inset on the right-hand side of (a) the position of the susceptibility maximum Tm vs D is depicted. 



z direction. Consequently, at zero temperature Xq~(0) = 2/D decreases with increasing D, 
and at intermediate temperatures Xq~{T) exhibits a maximum. 

Considering the longitudinal susceptibility Xo (Fig. Et>), it shows qualitatively the same 
behavior as in the D = limit; in particular, Xq z diverges as T — * indicating the ferromag- 
netic phase transition. For T >T the Curie law \t~ — ^Xo — 4/3T holds approximately, 
where, e.g., T = 2.5 (3) for D = 0.1 (1) and 5=1. 
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FIG. 7: Specific heat C for S = 1 vs T/D at D = 1,5,7.4, and 25, from right to left, comparing 
the Green-function (solid) with the ED (o, N = 12) results. The inset shows the positions of the 
maxima Tg in the temperature dependence of the specific heat as functions of D. 

C. Specific heat 

In Figs. [7| to the temperature dependence of the specific heat for the S = 1 chain is 
presented. As the main result, the ED on N = 12 chains with periodic boundary conditions 
yields two maxima at Tg and Tg , if D > Dq with Dq = 7.4, and only one maximum at 
Tg for D < Dq. Let us first consider the specific heat for D ^ 1 plotted in Fig. The 
Green-function results for N — > oo, agreeing with those for N = 12 within the accuracy of 
drawing, are in a very good agreement with the ED data. Our results for the maximum 
positions nearly agree with those of Blote^ obtained by the ED of N = 7 chains with open 
boundary conditions and subsequent extrapolations to N — > oo. For example, for D = 1 (5) 
we get Tg = 0.57 (0.85), as compared with Tg = 0.5 (1.0) in Ref. for D = 10 > Dq we 
obtain the maximum temperatures Tg ± = 0.85 and Tg 2 = 3.93 which are slightly larger than 
the values found in Ref. |l8|, Tg = 0.79 and Tg 2 = 3.77. At D = D the specific heat reveals 
a plateau within a small temperature region, 2.4 < T < 2.7 (cf. Fig. EJ). Correspondingly, 
the dependence on D of the maximum position exhibits a jump at D , as seen in the inset 
of Fig. For D > Dq, following the reasonings of Ref. Q the upper maximum at Tg 2 may 
be interpreted as Schottky anomaly due to the D term. 

To analyze the finite-size effects on the specific heat for D ^ 1 and the accuracy of the 
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FIG. 8: Specific heat for S = 1 at D = 1 and D = 25 (inset: N = 4, 12) and its dependence on 
the number of spins, where the Green-function results (line styles) are compared with the ED data 
(circles). 

Green-function approach in dependence on the chain length N and on D, in Fig. |H1 the 
specific-heat curves for different values of iV with 4 ^ iV ^ 12 and for D = 1 and D = 25 
(see inset) are plotted. As can be seen, the deviation of the Green-function results from 
the corresponding ED data decreases with increasing iV and D. Comparing the curves for 
D = 1 and D = 25, the finite-size effects decrease with increasing D. Moreover, they are 
found to decrease with increasing temperature which is not shown in Fig. |HJ where, e.g. for 
D = 25, only the low-temperature maximum is depicted (cf. Fig. Ej). 

Considering the specific heat at small anisotropies, the detailed analysis of our ED cal- 
culations for different chain lengths with periodic versus open boundary conditions reveals 
considerable finite-size effects, in contrast to the case D ^ 1 discussed above, and a re- 
markable dependence of the ED data on the chosen boundary condition. In this paper we 
prefer to use periodic boundary conditions, since, due to the translational symmetry imply- 
ing equivalent lattice sites, (i) the finite-size effects are expected to be less pronounced and 
(ii) ED calculations for larger systems (N ^ 12) can be performed, as compared with open 
boundary conditions used by Blote^ for N ^7. In Fig. El we illustrate the finite-size effects 
and the influence of boundary conditions for D = 0.1. The ED data for N = 4 yield a 
maximum at TJ? ~ 0.05 which vanishes for iV = 12. This may be understood as follows. In 



16 




FIG. 9: Specific heat for S = 1 at D = 0.1 and its dependence on the number of spins, where the 
ED results for N = 12 (o) and N = 4 (o) are compared with the Green- function theory for N — > oo 
(solid) and N = 4 (dashed). In the inset, the dependence of the specific heat on the periodic (solid) 
versus open (dashed) boundary conditions in the ED calculations for N = 8 is demonstrated. 

finite systems the spin excitations are gapped, even in the D = limit, where the finite-size 
gap Aat scales as N~ a with a > 0. If D < An, a low-temperature Schottky-type anomaly 
in the specific heat may appear and vanish for larger iV with D > An- Note that both ED 
curves approach each other at T > 5. In the Green- function theory for N = 4 a maximum 
is also found at the same temperature T£ = 0.1, but with a too large height. However, 
for N — > oo this maximum is only weakened, but does not disappear. In view of our ED 
results for N = 12, this behavior of the specific heat has to be considered as an artefact of 
the Green-function theory for small anisotropies. As can be seen from the inset of Fig. 
the use of open boundary conditions favors the appearance of a spurious low-temperature 
maximum in the specific heat. 

In view of our analysis described above, we consider the ED results by Blote^ on the 
specific heat of the 5 = 1 ferromagnetic chain with D < 0.25 as questionable, in particular, 
because the extrapolation of the data for small systems with N ^ 7 and open boundary 
conditions was performed. Our ED results on the specific heat at small anisotropies quali- 
tatively deviate from the data by Blote.— In Ref. Qtwo maxima were obtained not only for 
large values of D (see above), but also for D < 0.25, where for D = 0.1 a low-temperature 
maximum was found at = 0.12. In our ED data at large enough iV such a maximum 
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FIG. 10: Specific heat C of the Ni complexes NiBr2 • 2py (a) and NiBr2 • 2pz (b), where the Green- 
function theory (solid) is fit to the experimental data (♦, Ref. Il9r) . The insets show the predicted 
temperature dependences of the transverse magnetic susceptibility xti~ ■ 

does not appear. 

D. Comparison with experiments 

Finally, let us compare the results of the Green-function theory with some experiments on 
Ni complexes^ and derive predictions for quantities not yet measured. In Fig.^Jthe specific 
heat of the di-bromo Ni complexes NiBr 2 L 2 with L=pyrazole (pz, N2C3H4) and L=pyridine 
(py, NC5H5) is depicted. Those compounds can be considered as weakly antiferromagneti- 
cally coupled ferromagnetic chains with a large easy-axis single-ion anisotropy.^ The small 
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FIG. 11: Spin- wave spectra tOq for NiBr2 • 2py (solid) and NiBr2 • 2pz (dashed) at T = 5K 
and 10K, from top to bottom, and inverse correlation length (C 2 ^) -1 (inset), as predicted by the 
Green-function theory [q and £ zz are given in units of the lattice spacing.). 

values of the Neel temperatures T N indicated in Fig. [TUl reflect the pronounced quasi- ID 
behavior. The anomaly of the specific heat at T/v cannot be described by our theory for 
a purely ID system. For NiBr 2 • 2py (Fig. [TUk). this anomaly masks the low-temperature 
maximum at TJ? . At sufficiently high temperatures T > T/v the systems exhibit ID be- 
havior, and the theory may be compared with experiments. For NiBr2 • 2py (2pz) the fit 
to the specific heat data yields J = 0.4meV (0.48meV) and D = 3meV (2.7meV) so that 
D/J = 7.5 (5.6), where the first ratio slightly exceeds D /J. Note that those values nearly 
agree with the findings of Ref. 19. Using the fit values for J and D we calculate the tem- 
perature dependence of the transverse magnetic susceptibility xti~ = ^"b^aXo" {Na is 
the Avogadro constant). The results (see insets of Fig. [TUJ) show a maximum of Xm~(^) a ^ 



which should be confirmed experimentally. 

Furthermore, in Fig. ^2 we show the spin-wave spectrum and the correlation length 
(inset) calculated for the J and D values given above. Those results may be verified by 
neutron scattering experiments on single crystals. As disscussed in Sec. A, spin-waves in 
the paramagnetic phase may be observed, if q ^> {^ zz )^ 1 - For example, at T = 5K this 



T*>T N , where 



'-px 

± m 



5.35K ; 2py 
6.25K ; 2pz 



(37) 
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condition may be fullfilled for NiBr 2 • 2py (2pz) with (^ 2Z )" 1 = 0.47 (0.16). At T = 10K we 
have (£ 22 ) -1 = 1.16 (0.74) for the 2py (2pz) complex, so that only Brillouin-zone boundary 
magnons in NiBr 2 • 2pz may be observable. 



V. SUMMARY 



In this paper we have developed a Green-function theory for S ^ 1 ferromagnetic Heisen- 
berg chains with an easy-axis on-site anisotropy, where products of three spin operators 
are approximated in terms of one spin operator. Moreover, we have performed exact di- 
agonalizations of chains with up to N = 12 sites imposing periodic boundary conditions. 
To investigate the spin-wave picture in the paramagnetic phase, we have calculated the 
magnon spectrum and the correlation length. The thermodynamic properties (longitudinal 
and transverse susceptibilities, specific heat) at arbitrary temperatures were found to be in 
good agreement with the exact results for finite chains. A detailed analysis of the ED data for 
the specific heat yields two maxima in the temperature dependence for D/J> 7.4, whereas 
for D/J< 7.4 only one maximum appears. Our results at low ratios D/J contradict those 
of Ref. Q obtained on smaller chains with open boundary conditions. The Green-function 
theory was compared with specific heat experiments on di-bromo-pyrazole/pyridine Ni com- 
plexes, and predictions for the spin-wave spectrum, the correlation length, and the maximum 
in the temperature dependence of the transverse magnetic susceptibility were made. 
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